Equilibrium Dynamics of Microemulsion and Spon 

Phases 



G. Gompper and M. Hennes 
Sektion Physik der Ludwig-Maximilians-Universitat Miinchen 
Theresienstr. 37, 80333 Miinchen, Germany 

February 1, 2008 

Abstract 

The dynamic structure factor G(k, to) is studied in a time-dependent Ginzburg- 
Landau model for microemulsion and sponge phases in thermal equilibrium by field- 
theoretic perturbation methods. In bulk contrast, we find that for sufficiently small 
viscosity rj, the structure factor develops a peak at non-zero frequency lo, for fixed 
wavenumber k with /cq < k^q. Here, 27r/g is the typical domain size of oil- and 
water-regions in a microemulsion, and ko ~ rjq^. This implies that the intermediate 
scattering function, G(k, t), oscillates in time. We give a simple explanation, based 
on the Navier-Stokes equation, for these temporal oscillations by considering the 
flow through a tube of radius i? ~ vr/g, with a radius-dependent tension. 

PACS numbers: 61.20.Gy, 64.60.Ht, 82.70.-y 
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1 Introduction 



The behavior of self-assembhng amphiphihc systems has attracted much attention recently 
[|1], 0. The unique properties of these systems can easily traced back to the ability of the 
amphiphile to reduce the interfacial tension between two normally inmiscible fluids like 
oil and water by several orders of magnitude. One of the most interesting phases of 
ternary amphiphilic systems is the microemulsion, a homogeneous, isotropic mixture of 
all three components. It is now well established that microemulsions with medium- or 
long-chain amphiphiles consist of coherent oil- and water-regions, which are separated by 
amphiphilic monolayers. Similarly, the sponge phase in aqueous amphiphilic solutions 
consists of distinct water- regions, which are separated by an amphiphilic bilayer 
Thus, the static behavior of these systems is now relatively well understood. Much less 
effort has been made so far, however, to study the dynamical behavior of these systems. 
Most of these studies have considered non-equilibrium situations, typically spinodal phase 
separation ||, spontaneous emulsification [0, oil solubilization in aqueous surfactant 
solutions and amphiphile aggregation at an oil/water interface 0. 

We want to investigate here the dynamic structure factor G(k, cj), both with bulk 
and film contrast, of microemulsions in thermal equilibrium. It is well known from many 
studies near critical points that systems belonging to the same static universality class may 
show completely different dynamical behavior |TD, H]. It is thus essential to clarify which 
variables have to be included in a dynamic model of microemulsion and sponge phases. 
A microemulsion, for example, has two conserved densities, the amphiphile concentration 
^/'(r,t), and the difference of oil- and water-concentration $(r,t). A sponge phase, on the 
other hand, may also have two conserved densities, if the amphiphilic bilayers are almost 



inpenetrable to water, or only one, the amphiphile concentration, if they are not [12 



Nevertheless, the amphiphile concentration ip may not have to be considered explicitly, 
if its relaxation is fast compared to that of the order parameter $. This is the case we 
want to study in this paper. The opposite limit, where the amphiphile dynamics is slow 
compared to the order parameter dynamics, is more complicated and will be considered 
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elsewhere. We will show below that it is essential to include the hydro dynamic variables 
in a model for microemulsion and sponge phases. In general, these are the pressure p{r, t), 
and the longitudinal and transverse components of the momentum density, jL{r,t) and 
jT(r,t). It has been shown that in the vicinity of a critical point, p and Jl are irrelevant 
variables, and thus can be omitted in the model |T^. We will assume that the same is 
the case here. However, p and j^, are essential for studies of sound propagation in these 
systems |T3|, O, |15 . 



2 Ginzburg-Landau Model 



The time- dependent Ginzburg-Landau model we consider is a generalization of the Siggia, 
Halperin, and Hohenberg model H [[ 



T0[], consistent with linear hydrodynamics. Our 



hydrodynamic variables are the local order parameter $(r,t), which is proportional to 
the local concentration difference between oil and water, and the transverse component 
of the momentum density jT(r, t). In this case, the stochastic equations of motion can be 
written as 



'dt 
dt 



r$v2^ + (7oV(<i>^) + a 



5$ 



I 

5T\ 



where the Langevin forces Q have the usual correlations [IG 



< C<i.(r, t)C$(r', t') > = -2r$V'5(r - r')5{t - t') 
< CT,a(r, t)CT,/3(r', f) > = -2TTV^5iT - r')5(t - t')5, 

The equilibrium free energy functional has the form 

^= I d^r fc(V'$)2 + ^($)(v$)2-F/($) + ^j| 



0/3 



(1) 

(2) 



(3) 
(4) 



(5) 



It has been shown [17-22] that the static behavior of ternary amphiphilic systems is 
very well described by this free energy functional, with jT(r, t) integrated out, when the 



functions / and g are chosen to be 

gi^) = 6o + &2$' (6) 
/($) = ra*' + r4<l>' + rg*' , (7) 

with c > 0, 62 > 0, and rg > 0. Three-phase coexistence between oil-rich, water-rich and 
microemulsion phases is obtained for r4 < 0, r2 > and 60 < 0. As discussed above, we 
assume that the transport of the amphiphile is much faster than all the other transport 
processes. Thus, the deviation of the amphiphile concentration ip{v^t) from its mean 
value ip =< ip{v,t) > does not appear explicitly in our model. However, the average 
concentration ip and the properties of the amphiphile enter the model via the parameters 
of the functions / and g. The value of in the microemulsion, bo, for example, 

can be made negative by increasing the amphiphile concentration or the strength of the 
amphiphile [Q. Since we are interested here only in the behavior of the microemulsion 



phase, it is sufficient to retain only quadratic terms in the free energy functional (|^), 
so that 

g{^) = bo (8) 

/($) = ra*^ (9) 



The same model can also be used to describe sponge phases ||2J, |, In this case, 
the order parameter $(r,t) is identified with the local concentration difference between 
water on one side ("inside") and the other side ("outside") of the amphiphilic bilayer. 
Since we use a conserved order parameter in ([|) and (|), our results apply to sponge 
phases in which the leakage time of water through the bilayer is very large [|1^ . 

3 Van- Hove Theory 

When non-linear contributions in the Langevin equations (|l]), (0) are neglected, all struc- 
ture factors can be calculated easily [O. In this approximation, the Fourier transform of 
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Eq.(|I|) reads, for example, 

- ia;<l>(k, cu) = -2T^k'^ {ck^ + bok'^ + rs) $(k, cu) + C$(k, cu) . (10) 

Since (|10D is a linear relation between the random variable $(k, u) and the noise C$(k, uj), 
the noise correlation function (^) immediately implies that the dynamic water-water cor- 



relation function has the form 



a;2 + [T^k^xoik) 



^112 ' 



where 



Xo 



{k)'^ = 2 (cA:^ + bok^ + 



r2 



(12) 



is the static structure factor. For 6g < 4cr2, the static correlation function in real space 



is given by |T7 



where ^, with 



G$$(r) = — e"^''^ sin(gr) 
r 



^ 2V c 4 c 



(13) 
(14) 



is the correlation length, and 27r/g, with 

^ 2V c 4 c ^ ' 

is the characteristic domain size of coherent oil- and water-regions. It follows from Eqs. 



p!2D and (13) that there are three important lines in the phase diagram |17, 18]. With 



decreasing bo, the first is the "disorder line", 60 = V^^, where the correlation function 
changes its behavior from a monotonic decay to damped oscillations. The second line is 
the "Lifshitz line", bo = 0, where a peak in the structure factor begins to move to non-zero 
wavenumber. Finally, the line 60 = — -\/4cr^ is the spinodal to the lamellar phase. Note 
that the shape of the static structure factor is characterized by a single dimensionless 
parameter q^. 

The intermediate scattering function, ^^^(k, t), is obtained from Eq. ([Tl|) by Fourier 
transform, and reads ITT 



(16) 



Both Eqs. ( |TT| ) and (|TB|) can be used to define a characteristic frequency, 

uj, = r^k\oik)-' . (17) 

This characteristic frequency can be written in the scahng form 

LJc = k'no{k^,k/q) , (18) 

where the dynamic exponent z = 6 and the scahng function 

Qoix, y) = 2cT^ [l + 2{x-^ - y-^) + {x~^ + y-^f] (19) 

are obtained easily from Eq. (0). The asymptotic behavior of ujc in various hmits is 
hsted in Table I. The inverse scaling function, i.e. /uoc, is shown in Fig. 1. Two slow 
modes are present in a structured microemulsion, one for A; ~ 0, the other for A; ~ g. 
In the first case, the conservation of $ causes the divergence of the relaxation time in 
the limit — 0, as in any other binary fluid. In the latter case, critical slowing down 
occurs as the transition to the lamellar phase is approached; this slow mode is thus due 
to the increasing stability of fluctuations with the wavevector g, which characterizes the 
structure of the microemulsion. 



4 Perturbation Theory 



In order to treat the non-linear interactions, we employ a field-theoretic approach |^ 
using response fields $ and Jt- The two-point correlation and response functions 
G^-^^^{1,2) =< ipi{l)ipj{2) > are given exactly by Dyson's equation. 



G 



;i,2) 



:i,2) 



(20) 



in terms of the self-energies S^^,,/.^-, with ipi G {$, jr, jr}- Here, G^^^^^.{1, 2) denotes the 
correlation and response functions of the linearized theory. The inversion of the Dyson 
equation yields in particular the correlation function 

2r$fc^ + 



G#$(k, u) 



iu + r$A;2xo(fc) ^ - 



(21) 



In the one-loop approximation, the self-energy S^|,(k, cu) is given by 

^ s _ _ 2 /■ Xo(p) k ■ rp_k ■ k 

L,,(k, u)- g,J ^2^^3 ^^^^^ ^ ^^^p _ j^^^ ^ r^p2^^(p)-i ' 



with the projection operator 



The Feynman diagrams are shown in Fig. 2a. The second self-energy, T,^^, in Eq. ( piD 
can be obtained from Eq. ( p2D by the fluctuation dissipation theorem, which leads to 

E^^{k,u) = -2xo{k) Re{S^|,(k,cu)} , (24) 

or directly from the loop expansion, just as The calculation of the other correlation 
and response functions is described in the appendix. 



To calculate the characteristic frequency Uc in perturbation theory, we define p8 



_ G^^ik,t = 0) IZ.tG^^ik,u;) 2xo{k) 

dtG^^{k,t) G^^{k,uj = 0) G$$(fc,cu = 0) ■ ^ ^ 



In the case of the van- Hove theory, this definition agrees with Eq. ([T7|). 

The results of van-Hove and perturbation theory for Uc are compared in Fig. 3. The 
dependence of uJc on the wavevector k shown in Fig. 3a demonstrates that uJc becomes 
larger (i.e. the typical relaxation times becomes smaller) when the non-linear interactions 
are taken into account. This shows that the hydrodynamic fiow speeds up the relaxation 
process of the order-parameter correlations. This hydrodynamic relaxation mode must be 
a fiow of oil and water through their multiply connected networks which are characteristic 
for a microemulsion. 

Fig. 3b shows uJc for k = q. This demonstrates again the point we have made above 
for the van-Hove theory that there is a slow mode, with a relaxation time which diverges 
as the spinodal of the transition to the lamellar phase is approached, i.e. for l/{qC,) —>■ 0. 
However, it can also be seen from Fig. 3b that the relaxation time obtained from the 
perturbation theory remains finite in the limit qC, oo. This is a disturbing result, 
since it can easily be seen from Eqs. (^, (||) and (|^) that the static spinodal remains 
unchanged. 



It is not difficult to find the origin of tliis problem. It is the use of the van- Hove response 
function G'^?^ of the momentum density in Eq. ( ^9]) for the self-energy, compare Fig. 2a. 
In the van-Hove approximation, the response and correlation functions of the momentum 
density contain no information about the structure of the microemulsion, as can be seen 



from Eq. (|46|) . To avoid this problem, all one- loop correlation and reponse functions 
have to be evaluated self-consistently. Since this requires a considerable numerical effort, 
we only use the one-loop result for Gj^ and Gjj in Eq. (^Q]), which is the dominant 
contribution. With this modification, the characteristic freqency for k = q now correctly 
vanishes at the spinodal, as can be seen in Fig. 3c. 

The dynamic correlation function u) as a function of u for three different values 

of the wavevector k is shown in Fig. 4. Note that for A; ~ g/2 a peak develops for uj > 
in the limit of large g^, i.e. for strongly structured microemulsions. This peak indicates 
that the correlation function oscillates in time. Indeed, the Fourier transform of the data 
shown in Fig. 4b has temporal oscillations, see Fig. 5b. For k > q, the peak is much 
broader, so that there are still oscillations, but with a rapid decay, see Fig. 5a. Finally, 
for k q, the peak at finite u is absent, and the correlations function for large times 
decays monotonically, as can be seen in Fig. 5c. The peak at finite uj and the temporal 
oscillations of the dynamic correlation function are reminiscent of the peak at finite k 
and the spatial oscillations of the static correlation function. However, while the latter 



has been interpreted convincingly in terms of the microemulsion structure ||17| , P9|] , the 
mechanism of the former is unclear. 

Finally, we show the dynamic correlation function (^^^(k, cu) as a function of the 
wavevector k for fixed dimensionless freqency uj/T^ = 0.3 in Fig. 6. This function 
vanishes in the limit k ^ due to the conservation of the order parameter. 

5 Channel Flow 

To understand the temporal oscillations of the correlation function in microemulsions, we 
follow an idea which was used to explain the dynamics of spinodal phase separation in 
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the bicontinuous regime It is argued in that case that the relaxation is dominated 
by the flow through the channels; this leads to a growth law for the coarsening of the 
bicontinuous structure, which is linear in time t, rather than the t^^'^ behavior observed 
for droplets The force which drives the flow through the channels is the interfacial 
tension. 

To derive an equation which describes the time dependence observed in microemul- 
sions, we start from the linearized Navier-Stokes equation for the velocity field v(r, t), 

|-v = --Vp+^VV (26) 
ot p p 

where p is the mass density, and r] is the viscosity. Consider a spherical droplet of radius 
L within the microemulsion phase, in which the concentrations of oil and water deviate 
slightly from their average values, as shown schematically in Fig. 7. This implies that the 
cross-section of the oil- and water channels in this region is somewhat larger (or smaller) 
than on average. The easiest way for the system to get rid of the extra oil and/or water is 
to let it flow through the channels to the outside of the droplet. We now assume that the 
flow through this multiply connected network of a microemulsion is similar to the flow 
through a cylindrical tube of length L. For a cylinder in the z-direction, the flow field is 
approximated by Poiseuille flow, 

v(r,t) = a(t)[i?(t)2-r2]e, (27) 

where R{t) is the radius of the tube and a{t) describes the magnitude of the velocity. 
Further, we employ the Laplace equation to relate the pressure gradient to the tube 
geometry and to the interfacial tension cr, 

Vp ^ -^e, = --^e, (28) 

for a tube of length L. The Navier-Stokes equation ( PBD then implies that the average 
velocity (averaged over the cross-section of the cyhnder), 

1 i-m) 

v{t) = — - / v{r,t)dr (29) 



R{t) 
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is determined by 



d _ _R a 61] V 



where R = dR/dt. Note that it follows immediately from Eqs. (P7|) and (53) that 



V = ^aR . Since for an incompressible fluid 



we finally arrive at 



The analysis presented so far is not specific to bicontinuous microemulsions, but equally 
well applies to the case of bicontinuous phase separation. The specificity of amphiphilic 
systems comes in when we consider the interfacial tension a. It is well known that the 
interfacial tension between the oil-rich and the water-rich phase in these systems is very 



small. What we really need in Eqs. (p8|) and (|32D, however, is not the tension of the 
oil/water interface, but the driving force for a change in the tube radius. It has been 
shown in Refs. and |21| that in the lamellar phase the tension vanishes identically 



for the equilibrium separation of the interfaces. However, if the spacing between the 
lamellae is increased (decreased), a negative (positive) tension develops which drives 
the interfaces back to their equilibrium distance. We assume here that the same is true 
in the microemulsion. For small deviations from the equilibrium radius Rq, we can thus 
write 

a{R)=p{R-Ro) + 0{{R-Ro)^) (33) 
with a positive constant p, and 

R{t) =Ro + e{t) . (34) 
By inserting these results into Eq. (|32D, we find that e{t) is determined by 
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to leading order in e. Thus, e{t) satisfies the differential equation of a damped harmonic 
oscillator. The radius oscillates if the damping is small enough, i.e. if 

Thus, there should be oscillations if the viscosity r] of oil or water are small enough. We 
want to emphasize that the i?-dependence of the tension a in Eq. ( |55| ) is unrelated to 
a stretching of the amphiphilic monolayer (i.e. a change in the area per headgroup), 
since we are assuming in our model (^) that the amphiphile concentration samples its 
full equilibrium distribution on a time scale rapid compared to changes in the oil- and 
water-concentrations. 

With these results, we can now predict the behavior of the dynamic water-water 
correlation function. We identify the tube radius Rq with n/q, and the tube length L 



with vr/fc. The viscosity in our Ginzburg-Landau model is given by r/ = pTx [0- In 
terms of the variables q, k and Tt, the condition (|36|) then reads 

> aor|g^ , (37) 

where Qq = 18p/(7r^p). Here, both the parameter p and the mass density p depend on 
the coefficients c, bo and r2 in the free energy functional (|^), but not on any dynamic 
coefficients in Eqs. (|lD, (0). While the variation of p is probably small, p can vary con- 
siderably. Its functional dependence can be obtained from a calculation of the interfacial 
tension |^T|. In the following discussion, we will consider a fixed set of parameters of the 
free energy functional (^, (H), (H), so that Oq and q are constant. 

In case the inequality ([37| ) is satisfied, the solution of Eq. (pSf ) then implies that the 
asymptotic decay of the correlation function should have the form 

Gci>$(k, t) = B exp(-t/r) cos{iy{k)t + if) (38) 

where 



^^kf = '^e-a^Vlq'' (40) 

Co 
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with constants B, (p, and ai, 02, 03. We want to emphasize that the relaxation of regions 
with increased oil- or water-concentrations is facilitated by the flow mechanism on length 
scales larger than the diameter of a tube, i.e. for k <^ q. This is the regime where we 
expect our result ( p8| ) with (|37| ) to describe the dynamics correctly. On the other hand, 
on length scales smaller than the tube radius, i.e. for k > q, the dynamics should be 
determinded by undulations of the microscopic oil/water interfaces, or by the fluctuations 
within the channels, which are ignored in our derivation. 

Our predictions (|37D and (|38D with (|39D and (^OD for the correlation function can now 



be compared with the result obtained from the time-dependent Ginzburg-Landau the- 
ory. First, note that the inequality (|37| ) implies that there should be no oscillations for 
sufficiently small k, in agreement with the result shown in Fig. 5c. To extract the relax- 
ation time and the oscillation frequency from the perturbation theory, we fit the inverse 
correlation function, for small 00 (up to the peak position) to a fourth-order polynomial, 

to)-' = fioik) + fX2ik)tu^ + /i4(fc)cu^ + Oico'^) (41) 

and then use the equivalent of Eqs. ([I^ and (|15D to calculate r and z/. We have compared 
the results in a few cases with the explicit Fourier transform, G$$(k, t), and have found 
good agreement. The relaxation time r and the oscillation frequency u as functions of 
k and F^- are shown in Fig. 8 and 9. It can be seen that in the regime k <^ q there 
is semi- quantitative agreement with predictions ( P^ and (|D). The linear dependence 



of 1/r on Fy and of z/^ on F^ is found to be satisfied very nicely, see Figs. 8a and 9a. 
However, there is a considerable /c-dependence of r, see Fig. 8b, although from (|39| ) it 
is expected to be constant. In the regime k'^q, where we do not expect our arguments 



leading to Eqs. (|39D and (1^) to apply, the behavior of r and u as a. function k is indeed 
qualitatively different. The relaxation time is found to have a maximum at k = q, while 
the oscillation frequency almost vanishes, see Figs. 8b and 9b. Finally, for k > q, the 
relaxation time decreases rapidly. 
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6 Sponge Phases 

It has been pointed out in Refs. p4|, ^ that in sponge phases the order parameter <l> 
and its fluctuations cannot be observed directly, since water on one side of the membrane 
cannot be distinguished from water on the other side. The order parameter fluctuations 
can been seen indirectly, however, via the fluctuations of the amphiphile density ip{v,t). 

Although our model does not contain any explicit amphiphile degrees of freedom, 
we can nevertheless calculate an amphiphile correlation function, if we assume that most 
of the amphiphile is located at the $(r, t) = surfaces. Here, we assume again that the 
amphiphile concentration equilibrates instantaneously with relatively slow changes of the 
order-parameter field $(r, t). In this case, the amphiphile correlation function is given by 

Gfumi^, t) = No< <5($(r, t)) 5($(0, 0)) > , (42) 

with a normalization factor Nq, which is chosen such that limr^oo linit-*oo G^/i/m(r, t) = 1. 
The average can be calculated easily for Gaussian fluctuations, where one finds ^ 



No r 



-1/2 



Gfum{r,t) = ^ < $(0,0)2 ^2 _ < $(r,t)$(O,0)) ^" . (43) 

2,71 L J 

For not too small r or t, this result can be expanded to give 

GfuU^, t) - GfiiUr = oo, t = oo) ~ < $(r, t)$(0, 0)) . (44) 



Interestingly, in the static case Eq. ( p^ is just the result obtained from the two-order- 
parameter model, where the amphiphile density is included explicitly [^. We thus iden- 
tify the expression (Q) with the amphiphile correlation function, G^^{r,t). 

The results for the amphiphile correlation function, calculated with the van-Hove and 
the one-loop approximation for the order-parameter correlations, Eqs. (^Tj) and (P^D, are 



shown in Fig. 10 as a function of the frequency u for k = q/2. It can be seen that 
the scattering intensity in film contrast does not have a peak at finite u, but develops a 
shoulder at the same value of uj, where the order-parameter correlation function has its 
peak. 
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We want to point out that a very similar behavior is found when the expression 
is used to calculate the static film correlation function 0. Only when additional coupling 
terms in a free energy functional for the two order parameters $ and are taken into 
account, a peak at finite wavevector k appears [^. Thus, we expect that a time-dependent 
Ginzburg-Landau model for $ and ip will show a peak of ^^^(k, cu) at a finite frequency 

UJ. 

By a numerical Fourier transformation of the data of Fig. 10, we obtain the behavior 
of the intermediate scattering function in film contrast, G^^(k, t), as a function of time 
t. The result is shown in Fig. 11 for the same wavevector k = q/2 for which the order- 
parameter correlation function oscillates in time, compare Fig. 5b. We find again an 
oscillatory behavior, with a characteristic frequency which is roughly the same as the 
characteristic frequency u of G<i,$(k, t), but with a considerably smaller relaxation time. 

Note that formally there is a strong similarity with the spatial behavior of the static 
correlation, where one finds with Eqs. ([13D and (|^) 

G^4r)^^e-^'^/^sm{qr)\ (45) 

a function which shows the same characteristics as a function of r as G^^(k, t) does as 
a function of t. In particular, G^^(k,t) > G^^(k,t = oo). We want to emphasize that 
these results of the time-dependent Ginzburg-Landau theory agree with the arguments 
presented in Section 5. From these arguments, it follows immediately that G^^{k,t) 
should oscillate in the same intervall of wavevectors as the order-parameter correlation 
function itself. Also, the oscillation period of the two correlation functions should be the 
same. 



7 Summary and Conclusions 

We have studied the dynamic behavior of microemulsion and sponge phases in thermal 
equilibrium. For a time-dependent Ginzburg-Landau model, we have analysed the dy- 
namic structure factor G(k, u), both in bulk and in film contrast, which can be measured, 
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for example, by inelastic neutron scattering experiments in microemulsions. Our calcula- 
tion shows that for systems with a sufficiently small viscosity rj, or with a sufficiently large 
domain size 27r/q oi coherent oil- and water-domains, G$$(k, a;) for fixed wavenumbers k 
of order q/2 develops a peak or shoulder at finite uj. This peak imphes that the intermedi- 
ate scattering function, G$$(k, t), oscillates in time. To understand this surprising result, 
we have studied the flow through tubes in a very simple, cylindrical geometry, employing 
the linearized Navier-Stokes equation. Under the assumption that the interfacial tension is 
radius-dependent, and vanishes identically for the equilibrium radius of the tube, we were 
indeed able to reproduce the oscillatory behavior of (^^^(k, t). Furthermore, we found 
that the dependence of both the relaxation time and the oscillation frequency on the vis- 
cosity in the two approaches agrees very well, while the dependence on the wavenumber 
k is qualitatively correct. Thus, we believe that the existence of temporal oscillations in 
microemulsions is well established. 

The behavior of the amphiphile correlation function, G^^{k,uj), is found to be very 
similar. In the range of wavevectors, where G$$(k, o;) has a peak at finite frequency ujq, 
G^^CkjUj) develops a shoulder, which appears also at frequency cuq. The intermediate 
scattering intensity is found to oscillate in time. The characteristic time scale of these 
oscillations is the same in film and in bulk contrast. 

We want to emphasize that the two approaches we have used should apply in two 
different limits. The Ginzburg-Landau approach is based on the assumption that the 
order parameter fields vary slowly in space and time, i.e. it applies to the weak segregation 
limit. The hydrodynamic flow through a tube with hard walls, on the other hand, should 
describe the behavior of the strong segregation limit. Since the same behavior is found in 
both limits, it should also be found in intermediate cases. 

The mechanism, which leads to the oscillations of the intermediate scattering function 
in bulk and film contrast, can be summarized as follows. In a microemulsion or sponge 
phase, there is some prefered average distance between the amphiphilic membranes. If by 
thermal fluctuations this distance increases or decreases in some region of these phases, 
there is a force, the interfacial tension, which drives the system back to the thermal 
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average. It is negative (positive) for large (small) distances between membranes, and thus 
tends to increase (decrease) the interfacial area. A restoring force itself is not sufficient to 
produce oscillations. We also need a moment of inertia, which drives the system through 
the point where the force vanishes; in our model, this term is provided by hydrodynamic 
momentum conservation. And we need a sufficiently small friction coefficient; in the 
microemulsion or sponge phase, this can be achieved by a sufficiently small viscosity, or 
by a sufficiently large diameter of the tubes, through which the flow moves back and forth. 
Thus, oscillations can only be observed for strongly structured and swollen microemulsion 
or sponge phases. 

We want to conclude with a short dicussion of the dynamic behavior of the amphiphile. 
We have already pointed out that the radius- dependence of the interfacial tension is not 
due to a stretching of the amphiphile film, since in our model the amphiphile always 
samples its equilibrium configurations. Let us now consider the opposite limit where the 
dynamics of the amphiphile is very slow. When the oil- or water-concentration in a drop 
of radius L ^ Ti/q decreases, the interfacial area must increase. Since the amphiphile 
dynamics is slow, its concentration in the drop can be assumed to be approximately con- 
stant. This implies that the amphiphile density within the monolayer decreases, leading to 
an increase of the area per headgroup and thus to a positive interfacial tension. However, 
this is the situation where we already have a positive tension in our model (|T]), (^, (^ 
with an equilibrated amphiphile. Thus, if the amphiphile dynamics is slow, the stretching 
of the monolayer enhances the radius dependence of the interfacial tension, and thereby 
the oscillations. 
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8 Appendix: Correlation and Response Functions 



The inversion of the Dyson equation leads to the correlations functions cj) given 

in Eq. (|21|) and 

Gni^. ^) = , r; -^ „ , (46) 

as well as to the response functions 



and 



G'|,^(k,u;) 



-iuj + r$Pxo(A;) ^ - s^i, 



1 



(47) 



(4J 



The correlation and resonse functions, G^^^^^. (k, uj) of the linearized theory can be obtained 
easily from these expressions by ignoring the self-energy contributions. 
In the one-loop approximation, the self-energy T,^^(k,u) is given by 



E^l,(k,t^) 



9o / TTT-^ / t;- k ■ Tp_k ■ k 



(27r)3y 27r 



~Gf^{p-k,u-u')Gfl{p,u')^ 



see Fig. 2a. The explicit form has already been presented in Eq. 
The one-loop expression for S^j is given by 



S,j(k,..) 



2 / d^p f duj' 



(27r)3 7 27r 



Xo(p)-' - Xo(p - k)-i 4°i(p,a;')G'i"i(p -k,uj- J) 



(49) 



(50) 



see Fig. 2b, or equivalently [|T^ 



%(k,c.) 



-^0 



[xoip) ^ - XoiP - k) ^] Xo(p - k) 

P ■ ■ P 

-iu; + T^p'^xoip)'^ + r$(p - k)2xo(p - k)-i 



(51) 
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Finally, the self-energy, Ej--, in Eq. can again be obtained from Eq. (^Tj) by the 

fluctuation dissipation theorem, which leads to 



E-(k,^) = -2Re{s^.j(k, 



(52) 



or directly from the loop expansion, just as T,---. 
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regime 




7 J- 1 






2cr<,(r' + g')'A;2 




2cr$e"^A;2 


< A; < g 


2cr$54A;2 


= 0, A; ~ g 


8cr$?4(A; - qf 



Table /: Asymptotic behavior of the characteristic frequency uUc in the van-Hove approx- 
imation. 



22 



Figure Captions 



Figure 1: The inverse scaling function Q^lk^, k/q) ^ = /ujf. in the van-Hove approx- 
imation. 

Figure 2: Feynman diagrams which contribute to one- loop order to the self-energies 
(a) and (b) T,---. Here, the straight lines represent the physical fields, the wavy 
lines the response fields. The solid lines denote order parameter propagators, the lines 
with bars momentum propagators. The diagrams contain the vertex contributions u = 
PXo(p)"^ - kxo(k)"^ and V = p [xo(p - k)~^ - Xo(p)"^]- 

Figure 3: The results of van- Hove (dashed lines) and perturbation theory (full lines) 
for uOci (a) as a function oi k/q for fixed = 8.888, (b)(c) as a function of q^ for fixed 
k = q. (c) The behavior of Uc for large q^ in detail. Here, the dashed-dotted line is 
calculated with the one-loop approximation for Gjj and Gjj. The parameters in Eqs. (|1|), 
(I), (I), i) and (D are = 1, Ft = 1, = 2, c = 1, and = 1. 

Figure 4: The dynamic correlation function (^^^(k, cu) as a function of the scaled 
frequency cu/F^ for three values of the wavevector, (a) k = 1.3q, (b) k = q/2, and (c) 
k = q/10. The dashed line is the van-Hove result, the full line the result of the perturbation 
theory. The parameters in Eqs. ([|), (||), (|]), (§) and (|]) are F$ = 1, Ft = 0.3, go = 5, 
c = 1, r2 = 1, and q^ = 8.888. 

Figure 5: The intermediate scattering function (^^^(k, t) as a function of the scaled 
time tF$ for three values of the wavevector, (a) k = l.3q, (b) k = q/2, and (c) k = q/10. 
The parameters are the same as in Fig. 4. 

Figure 6: The dynamic correlation function G^^ik^u) as a function of the scaled 
wavevector k/q for uj/Tq, = 0.3. The dashed line is the van-Hove result, the full line the 
result of the perturbation theory. The parameters are the same as in Fig. 4. 
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Figure 7: A spherical droplet of radius L (marked by the dashed line) within the 
microemulsion phase, in which the concentrations of oil (dark grey regions) and water 
(white regions) deviate from their average values (schematically). 

Figure 8: The relaxation time r of the asymptotic form ( pSj ) of the intermediate 
scattering function (^^^(k, t), as obtained from pertubation theory, (a) 1/r as a function 
of Tt for k = q/2 (full line) and k = q/?) (dashed line), (b) r as a function of k for 
Ft = 0.2. The parameters in Eqs. ([l|), (|), (|), (|) and (D are F$ = 1, Ft = 0.3, go = 5, 
c = 1, r2 = 1, and qC, = 19.97. 

Figure 9: The oscillation frequency u of the asymptotic form ( ^8| ) of the intermediate 
scattering function G$$(k, t), as obtained from pertubation theory, (a) z/^ as a function 
of Fy for k = q/2 (full line) and for k = q/3 (dashed line), (b) z/^ as a function of /c^ for 
Ft = 0.2. The other parameters are the same as in Fig. 8. 

Figure 10: The dynamic amphiphile correlation function (^^^(k, u;) as a function of 
the scaled frequency tu/F^, for wavevector k = q/2. The dashed line is the van- Hove 
result, the full line the result of the perturbation theory. The parameters are the same as 
in Fig. 4. 

Figure 11: The intermediate scattering function in film contrast, ^^^(k, t), as a func- 
tion of the scaled time tF$, for wavevector k = q/2. (a) van- Hove result (dashed line) 
and the result of the perturbation theory (full line), (b) A magnification of G^^(k,t) 
(perturbation theory) on intermediate time scales. The parameters are the same as in 
Fig. 4. 
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